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Abstract 

Using the J\-Ji model, we present a description of quantum phase transition 
from Neel ordered to the spin-liquid state based on the modified spin wave 
theory. The general expression for the gap in the spectrum in the spin-liquid 
phase is given. 
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There has been recently considerable interest in magnetically disordered states in quan- 
tum spin models. Much of this interest stems from the connection of this problem to high-T c 
superconductivity. The ground state of undoped compound has long range antiferromag- 
netic order. It is well described by the Heisenberg model and has been studied by numerous 
methods.0 However introducing a small number of holes leads to destruction of long range 
order. The resulting state is still not fully understood. 

Destruction of long range order can be studied by introducing some frustration into the 
Heisenberg model. We will focus on the simplest possible model of such kind which is the 
J1-J2 model defined by 

H = Jx]TS r -S r , + J 2 S r-Sr'. (1) 

NN NNN 

In this Hamiltonian, the J\ term describes the usual Heisenberg interaction of nearest neigh- 
bour spins [S = |) on a square lattice, while the J 2 term introduces a frustrating interaction 
between next nearest neighbour sites. The J1-J2 model itself is hardly applicable to any 
real materials! (although it was originally proposed! to describe high-T c superconductors). 
However, it is valuable for demonstrating how long range order can be destroyed. For con- 
venience, we set Ji = 1 and denote a = J2/J1 (the notation of Ref. [I] is used whenever 
possible). 

For small a, the ground state is Neel ordered. For large a the system is decomposed 
into two Neel ordered sublattices which, however, have the same quantization axis. This 
is the so-called collinear state. Whether or not the Neel and collinear states are separated 
in parameter space by a state without long range order has been the subject of many 
discussions. Besides the many spin wave calculations,!! - ! the model has been studiedby 
the Schwinger boson mean field theorvB analysis of small lattices," a series expansion,^ a 
mean field theory of bond operators,@ and other methods!! Despite the numerous efforts, 
a strict answer has not been obtained. Nevertheless, since only the mean field spin wave 
theory and the essentially equivalent Schwinger boson mean field theory predict a first order 
transition from the Neel to the collinear state, while all other methods provide support for 
the existence of a different intermediate state, the latter scenario appears far more probable 
(this is also supported by a recent calculation of corrections to the mean field solution!) . In 
this work we assume that the system undergoes a second order quantum phase transition 
at a certain a = a c from the Neel to a spin-liquid state. 

For the Neel phase we will use the spin wave theory which gives the following description. 
The staggered magnetization = | (0 1 |0) | is equal to 
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- £sinh 2 tf k , (2) 
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where N is the number of sites on the lattice, the summation is performed over the Bril- 
louin zone of one sub lattice (\k x \ + \k y \ < 7r), and 6* k is the parameter of the Bogolubov 
transformation determined in the linear spin-wave theory^ (LSWT) by 

tanh2# k = -— -, (3) 

1 + a{r] k - 1) 

or in the mean field spin wave theory theory! (MFSWT) by the selfconsistent solution of 
Eq. (H) together with 
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" (mt + ^i) + a(mt + (7a)(%-l)' 
2 1 2 
Si = jyE 2 7k Sinh 2 ^ k ' ^ 2 = N ^ ^ k Sinh2 
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We have defined 



7k = - (cos k x + cos k y ) and rfc = cos k x cos fcj,. 

Numerical values of rrv in LSWT, MFSWT and other approximations are shown in Fig. [3]. 
Further, the dispersion of the Goldstone spin-wave excitations is 



2 (1 -tanh 2 2fl k ) 1/2 , LSWT; 
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'I - tanh 2 29 u ) 1/2 , MFSWT. ^ 



for k <C 1 is ek = c/c. We calculate Z c = c/cq (co = a/2 is the spin wave velocity at a = 
calculated by LSWT) 

J (1 -2a) 1 / 2 , LSWT; 
Zc = \ 2v/2(mt + gy) (| - a^±j) V ' , MFSWT. (5) 

The numerical value of Z c is given in Fig. 0. Note the remarkable agreement of the mean 
field solution and the 1/S expansion when the latter is converging. 

Let us now examine the accuracy of the spin wave theory as the staggered magnetization 
m) decreases. The spin wave theory starts with transforming the spin operators S r to 
bosonic operators a[ and a r using the Dyson-Maleev or Holstein-Primakoff transformation. 
In either case, S* = \ — n r with h r = a\a r (a spin up sublattice is considered). The physical 
states are those with n r = and 1. However, the Hamiltonian is in fact simplified so that it 
connects the physical states to unphysical (see discussion in Ref. |1|). To estimate the amount 
of unphysical states introduced to the wave function we calculated averages of higher powers 
of the n r operator 

(0|n r |0) = m = - - m\ (0|n 2 |0) = m + 2m 2 , (0\n 3 r \0) = m + 6m 2 + 6m 3 . (6) 

where m is defined by Eq. (Q). If only n r = and 1 were present, we would have n\ = n r 
(/ > 1). The problem of unphysical states becomes serious as m increases. To estimate the 
weight of the states with n r > 1, we proceed as follows. Expand the state obtained in the 
spin wave theory in states with a definite number of bosonic excitations at site r 

oo 

|0> = J» r >, (7) 

n=0 

where n r \n r ) = n r |n r ). Obviously, the correlators of Eqs. (||) can now be expressed as 



oo 

(0|f4|0> = $>'c 2 , 1 = 0, 1, 2, 3, 



(8) 



TABLE I. The weight of unphysical states when the series (0) is truncated at C4. 
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a The m for a = 

Now we truncate the series (0) at c 4 and with the truncated series we solve the first four of 
Eqs. (§) taking the left-hand-side from Eq. (|6]) . The results are presented in Table |. Because 
a truncated series is used, negative weights sometimes appear. However, for the purpose of 
estimating, the method is adequate. We see that at m — 0.5 unphysical states constitute 
about 20 % of the wave function and we conclude that the spin wave approximation should 
be quite reasonable even at this point. (However, we must restrict ourselves to considering 
only low powers of the operators aj and a r since high powers have large contributions from 
unphysical states and may lead to absolutely unphysical results.) 

We assume that the system undergoes a second order transition into a liquid state at a = 
a c (the point when the sublattice magnetization becomes zero). In terms of the initial (a = 
0) ground state, zero sublattice magnetization means that the ground state is a condenstate 
of many spin waves ak and 6k- To describe the emerging phase we must take into account 
their non-linear interaction. We cannot do this exactly. However, there is an approximate 
method using the suggestion made by Takahashi0 for the Heisenberg model at non-zero 
temperature. 

Following Takahashi, we impose an additional condition that sublattice magnetization is 
zero 



<0|S1 - S* B \0) = (0|| - ala r + \ - 6j6 r |0) = 0, 



(9) 



where A and B are the spin up and down sublattices. In essence it means an effective cut-off 
of unphysical states. In fact, Eq. @) together with conservation of the z-component of the 
total spin is equivalent to 
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(0|E4a k |0) = (0|E^k|0) = -. 



(10) 



The total number of single spin wave states in the magnetic Brillouin zone is N/2. Therefore, 
after introducing Eq. @ the effective number of allowed states in the Hilbert space of the 
system is 



( (N/2)\ 
\(N/4)\(N/4)\, 
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so that with logarithmic accuracy the correct dimensionality {2 N ) is restored. 

The constraint @ is introduced into the Hamiltonian via a Lagrange multiplier A. Now 
we we must diagonalize 
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H X = H-X E4« r + E^ • (11) 



B 



The simple (linear) second term here [taken together with Eq. (|9])] takes account of non- 



linear interaction of spin waves. Diagonalizing Eq. (|TT|) we get the spectrum of excitations 



which has a gap A oc v^A- When A 1, the spectrum is 

£ k = (A 2 + 4) 1/2 , (12) 

where e k is the dispersion from Eq. (|4]) and A is determined from (for clarity, we present 
now equations only in the form they have in LSWT) 

- , + 2ff 7 5,fi-iU0 (13) 



(2ir) 2 Ve k B k 

[here m) < is calculated from Eq. (0)]. The integral in Eq. flT3| ) converges at small k where 
e k « cfc. After integration we have for a > a c 

A = 7 r|m t |c 2 . (14) 

The correlation length in this state is £ oc 1/A. 

In the vicinity of the transition A = 7rc 2 A(a—a c ), where c c is the spin wave velocity at the 
critical point (c c m 0.71 whether LSWT or MFSWT is used, see Fig. 2) and A = —dm* /da 
is the slope of dependance of m* versus a (Fig. 1). We can estimate A from LSWT or the 
improved MFSWli and we obtain A m 3.3(a — a r ). 

It is easy to generalize the consideration to non-zero temperatures. In this case the gap 
A will be determined from 



m ] + 2 
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(2vr)He k E k — \2TJ\ °' (15) 

Similarly to the T = case, the integration can be performed exactly for T, A 1 yielding 

T ( A \ 
■log(2sinh — ) = -p„ (16) 

where 



2vr ° V 2T 



1 t 2 
H 4 

Equation (^) has been obtained in the Na — oo limit of the non-linear a-mo de]@ (iV d = 3 
is the number of components of the order parameter) which describes long range behaviour 
of spin systems using several nhenomenological parameters. Solutions of Eq. (|i~6"D in different 
regimes have been discussedcJ (see also Ref. |T^) and we will not dwell on this issue. 

To summarize, we have presented a description of quantum melting of long range an- 
tiferromagnetic order in the frustrated Heisenberg model. The suggested approach can be 
used for many systems, in particular to the doped t-J model.0 

We would like to thank F. Mila for drawing attention to the topic. We are also grateful 
to C. Hamer, M. Yu. Kuchiev, D. Poilblanc, and T. Ziman for valuable discussions. This 
work forms part of a project supported by a grant of the Australian Research Council. 
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FIG. 1. Staggered magnetization m) as a function of a in the Neel and collinear states. Dotted 
line is the LSWT result,! solid line the MFSWT result,i long dashed line the improved MFSWT 
result,! dashed line obtained by the 1/8 expansion! in the first (o) and second (•) order. 



FIG. 2. Spin wave velocity renormalization Z c as a function of a. Dotted line is the LSWT 
result,! solid line the MFSWT result,! dashed line obtained by the 1/S expansion!! in the first (o) 
and second (•) order. 
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